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^ : Abstract 

> 

The free energy costs for various defects within an Abrikosov lattice of 
vortices are calculated using the lowest Landau level approximation (LLL). 
Defect solutions with boundary conditions for lines to meet at a point (cross- 
CTS . ing defect) and for lines to twist around each other (braid defect) are sought 

for 2, 3, 6, and 12 lines. Many results have been unexpected, including the 
nonexistence of a stable two- or three-line braid. This, and the high energy 
cost found for a six-line braid lead us to propose that the equilibrium vortex 
state is not entangled below the irreversibility line of the high- Tc supercon- 



O . ductors or in a large part of the vortex-liquid phase above this line. Also, the 

' solution for an infinite straight screw dislocation is found, and used to give a 

limiting form for the free energy cost of very large braids. This depends on 
^ ' the area enclosed by the braid as well as its perimeter length. 

.<^ ; PACS no.s: 74.60.Ge; 74.20.De 
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I. INTRODUCTION 

There has been a lot of interest in the effects of thermal fluctuations on the behaviour 
of flux-lines ever since the discovery of the high-Tc superconductors. The high tempera- 
tures and fields at which these new materials remain in the mixed state, and their small 
coherence lengths, mean that these fluctuations may be strong enough to totally alter the 
nature of the vortex state from the mean-field solution of an Abrikosov lattice. Debate and 
controversy has continued on the possibilities of different vortex phases above and below the 
irreversibility line (at which the vortices become depinned). We believe that the results of 
the calculations in this paper shed some light on certain questions on the vortex state. Our 
primary conclusion is that there is very little entanglement near and below the irreversibility 
line. 

In this work, we start with the triangular vortex-lattice, and consider the free energy costs 
of topological defects within the lattice. Because of thermal fluctuations, the spontaneous 
formation of defects of flnite free energy cost is important in the description of the vortex 
state. These defects may include 'crossings' where two or more lines may pass through the 
same point, and 'braids' where a collection of lines twist around each other as one moves 
along the fleld direction. Our approach is always to begin with the Abrikosov lattice of 
vortices that corresponds to the mean-fleld solution for the Ginzburg-Landau free energy 
functional within the lowest Landau level approximationtla. We realize of course that the 
LLL approximation only holds strictly in the high fleld regime of the superconducting phase 
diagram, just below Hc2i and even that the existence of this lattice may be questioned. 
Having set up the perfect lattice of vortices, we than allow a restricted number of the 
vortex-lines to move subject to certain boundary conditions that deflne a topological defect, 
and then flnd the minimum free energy cost of such a defect. 

The free energy cost of a crossing conflguration provides an estimate of the energy barrier 
for vortex lines to cross through each other, or to cut and reconnect. Accurate estimates 
of the crossing energy barrier in the different regions of the H-T phase diagram are desired 
because of their relevance to the vortex dynamics of high-Tc superconductors. The energy 
barrier, [/x , to crossing/cutting of vortex lines is a fundamental parameter of the entangled 
stateScl. This has encouraged calculations of this energy in different approximationa§Q. The 
calculation in this paper is the flrst to consider the effects on Ux of the surrounding triangular 
lattice. 

The braid defects are examples of the sort of conflgurations one expects to flnd in the 
entangled vortex state. As well as affecting the response to pinning centers, they will have 
interactions with any component of current parallel to the applied fleld. The growth of 
braid defects may be a source of dissipation of such longitudinal currents. An important 
result of our calculation is that there are no stable two-line or three-line braid defects, at 
least within our approximations. The smallest braid defect that is stable, but which has 
a relatively high free energy cost, is the six-line braid. The non-existence of stable braid 
conflgurations with small free energies has the consequence that in a sample of typical size 
of high-Tc superconductor there is a large region of the H-T phase diagram at which there 
will be no braid defects present in thermal equilibrium. We estimate that there will be no 
braids for values of the reduced temperature far above the irreversibility line and quite close 



to the Hc2 line (see Section |VA|) . If there are no braid defects present in a given sample then 



the vortex system is disentangled. This is contrary to the picture of an entangled vortex 
hquid suggested by NelsonQ. 

In Section || we describe the formulation of our calculations. Section ^T^ contains our 
calculations and results for crossing and braid defects involving two, three, six and twelve 
lines. In Section |IV| we extend our method to large scale defects of the Abrikosov lattice. The 
free energy cost of an infinite straight screw dislocation is calculated, and the result is used 
to find a limiting form for the free energy cost of very large braids. This free energy is found 
to depend on the area enclosed by the braid as well as the perimeter length. In Section it 
is shown how our results for small braids lead to the conclusion that entanglements are not 
important in a large region above and below the irreversibility line of high-T^ superconductors 
of typical sizes. Other applications of the calculations are also suggested in Section 0. 

II. METHOD OF THE CALCULATIONS 

For the sake of clarity, and to set the notation used, the general method behind the cal- 
culations is now described. First, the Ginzburg-Landau theory is formulated and the Lowest 
Landau Level is defined. Then the ground state corresponding to an Abrikosov triangular 
lattice is presented, and it is shown how to construct deviations from the perfect lattice 
which remain in the LLL. The thermodynamics of fiux-lines within a type-II superconduc- 
tor are well described for fields H near Hc2{T) by Ginzburg-Landau theory. Within this 
theory, the free energy density is expanded in terms of the superconducting order parameter 
and its derivatives, 
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where b = V x A is the microscopic flux density. This equation applies for general anisotropy 
but only applies to a homogenous system (layering effects are ignored). For relevance to 
layered superconductors, we consider a limited anisotropy with nix = '^y = ^T^ab and niz = 
rUc- For a uniform external field Ho parallel to the z-axis, it is the Gibbs free energy, 
Q{ijj} = Fq + J-'{ip} — b.Ho that controls the properties of the systema 

Our main approximation (valid for high fields) is to restrict the order parameter to the 
'Lowest Landau Level' (LLL) subspace, defined by 



n^ = o. 
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where 11 = 11^,. -|- illy, U^ = —ih{d/dx) — 26^4^,., and 11^^ = —ih{d/dy) — 2eAy. Within this 
restriction, substitution into the Gibbs free energy leads to an effective free energy density 
oil 
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where au = a + {eh/mab)fJ'oHQ, and Pk = P — f'oi^-h/mab)'^- Here, an is our reduced 
temperature variable, which is negative below the Hc2 line and positive above. It is useful 
to make the transformations 




if) ^ IP = 

z^h= (^nzpiinj ^. (2.4) 

Thus, h is the distance along the c-axis in units of the mean field correlation length in this 



direction, ^c = y ^ /2mc|a_H-|. Substituting (|2.4|) into ( p.3|) and integrating over all space 
(for an infinite bulk superconductor) gives a total free energy 
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A consequence of using the LLL approximation is that solutions of ( |2.2| ) have the general 
form in the plane perpendicular to the magnetic fieldcl 

V^o(x,y) = /(^)e-'^^'/^ (2.6) 

where z is not the third dimension in space, but the complex variable, z = x + iy, and 
1] = $o/-B, {B = (b) is the mean magnetic fiux density, and $o = h/2e is the fiux quantum). 
f{z) is an analytic function of z. Any function of the above form is a function within the 
LLL subspace. f{z) may be written quite generally in a product form 

f{z)^l[iz-z,), (2.7) 

i 

where the zeros in f{z), z = Zi, correspond to the positions of vortices within the supercon- 
ductor. Of course it is well known that the order parameter that minimizes the free energy 
( p.3| ) in the mixed state is a triangular periodic lattice of fiux-linesa. The corresponding 
function f(z) for such a lattice is a Jacobi theta functionlll 

M^, y) = 0(r|O) = Ce-^t^s (y , ^) . (2.8) 

Properties of 0(r|O) and related functions 0(r|ro) that span the LLL subspace are given inS. 
The function |0(r|O)| has the periodicity r/ and vn where r/ = /(1,0), vn = ^(1/2, \/3/2), 
and r = xn + Wii = 1/2 + il\/?)/2. I is the distance between neighboring zeros (vortices) 
in (/)(r|ro). The area of the unit cell must contain one quantum of fiux, so that r] = Z^v^/2. 
The prefactor in ( |2.8| ) is found to be C = 3^^^, if we use the normalization conditionlil 

(1^01')=/ ^|0(r|O)|^ = L (2.9) 

J unit cell Tj 

To find the free energy of the Abrikosov lattice, we substitute ijj{x,y,h) = KtljQ{x,y) into 
( p.5|) . As l^/'ol is periodic, we can integrate over one unit cell to find the average free energy 
density 

/=(il)'^/ ^{-A-|,,„|VA-W,r}. (2,10) 

V 2mr I Pk J unit cell Tl '- -• 



This is minimized with respect to K (using (|2.9|) ) by the condition K"^ = 1/(3a, where (3a is 
the Abrikosov parameter (3^ = (I'^ol )/(l'^o| ) — 1 ■ 1596. The free energy per flux-hne per 
unit length for the Abrikosov lattice is therefore given by f^ = —an"^ /2[3k[3a- 

As any LLL function has the general form (p.7|), including f{z) = ^3 (vrz//, vrr//) (see 
Appendix ^, we can therefore construct a function representing the order parameter for a 
triangular lattice with one flux line displaced from z = zq to z = (q, hj 

(z — C ) 
i^ dispiacedir) = ipoir)- ^. (2.11) 

[Z — Zo) 

This cancels out the first order zero in ipQ at z = zq and replaces it with a first order zero at 
z = Co- Because this new function has the form of ( |2.7| ), it is still within the LLL subspace, 
even though it is obviously not the configuration of lowest free energy (which is ipo). (A 
similar method to this was used by Brandttj). (|2.11|) can be generalized to any number of 



displaced lines, to form order parameters for various defects of the perfect lattice within the 
LLL subspace. The free energy cost of any such defect is given by integrating over all space 
the difference in the free energy densities of the defect and the ground state: 
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In all the problems described in this paper, the procedure is to define a defect in the 
perfect lattice by allowing the positions of some of the vortices to vary with h subject to 
certain boundary conditions leaving the rest of the vortices fixed in their positions in the 
triangular lattice. The total free energy of the defect is then minimized with respect to 
the positions of the chosen vortices, subject to these boundary conditions. This is done by 
expanding the free energy density as a polynomial in a set of variables that describe the 
coordinates in the x-y plane of the chosen lines, say {Ci{h)} for the vortices labelled by i. The 
coefficients in the expansion will be functions of x and y only (as there is no h dependence 
of ipo). These coefficients must be integrated over the x-y plane to give the free energy per 
unit length along h, that depends on {Ci{h)} and their first derivatives. The form of the 
variables {Ci{h)} that minimizes the total free energy will be solutions to Euler-Lagrange 
non-linear differential equations. By solving these the lowest free energy cost of a given type 
of defect can be found. 

One of our main approximations is that we do not allow the lines not directly involved 
in the defect to move in response to the presence of the defect (i.e. the relaxation of the 
surrounding lattice). This may seem like a poor approximation that will seriously over 
estimate the free energies of the real defects. However, we have found that this is not the 
case when we allow the nearest neighboring lines to move (as in Section 1IIB| ) and we believe 



that relaxation of the lattice will only slightly reduce our values of defect free energy costs. 
The reason for this is that we are considering a lattice of vortex lines, which cost energy to 
tilt with respect to the field direction, so if we consider a localized defect taking place over 



a small length scale, L, then the surrounding lines would have to tilt considerably over the 
distance L if they are to reduce their interaction energy with the defect lines, and this will 
cost too much tilt energy. 

An important general result of our method is that for any defect of the ground state 
made from changing the positions of n of the vortices, the free energy cost will diverge 
logarithmically with the total size of the system, unless the n lines move symmetrically 
about their 'mean midpointtiy. That is, if n lines move from their ground state values {zi] 
to the new positions {Ci{h)}i then the condition for a non-divergent free energy is: 

n n 

Y.ah)=Y.^^^ (2-14) 

i=l i=l 

i.e. if we define the 'center' of the defect as the vector sum of the coordinates of the ground 
state lattice position, z,, then this condition says that the vectors of the coordinates of the 
defect line must always sum to this 'center'. An important point to note is the restriction 
in this derivation to a defect of n lines. If one wanted to consider a defect that by definition 
could not satisfy the symmetry relation ( |2.14|) then a finite free energy could still be obtained. 



but only by allowing the 'relaxation' of other lines in the system. The motion of these extra 
lines would be so as to ensure the symmetry condition for the total number of moving 
lines, and through this cancel out the logarithmic divergence. However, if it is possible for 
a simple n-line defect to satisfy ( p. 14 ), then the free energy cost will depend only on the 



locally surrounding lattice. We take advantage of this result in the following calculations to 
reduce the number of free parameters used for describing a given defect. 

III. CALCULATIONS FOR SMALL DEFECTS 
A. Two Lines Crossing 

The first defect to which the above method is applied is the case of two neighboring 
lines within the Abrikosov lattice moving together to meet at a point. The importance 
of this defect is that it is the configuration providing an energy barrier to the 'cutting' 
or 'reconnecting' of two vortices. The simplest way to estimate the energy of two lines 
crossing is to keep all the surrounding lines in their original positions in the lattice, and to 
assume that the two lines move symmetrically towards each other (see Fig [ip. This avoids 
a divergence in the free energy cost, and results in an order parameter depending on one 
parameter only. 

The calculation for this simple defect is described in more detail than the other problems 
in order to demonstrate the general procedure of these calculations. If we take the origin, O, 
of the complex z-plane to be at the midpoint of the two lines, then we have two vortices (= 
zeros in the order parameter) which move from z = ±1/2 to z = ±a{h), with the boundary 
conditions a(±oo) = ±1/2. Using ( p.llj ), the LLL order parameter for two lines crossing 
may be written 

Az ± a(h))(z — a(h)) , , 

^2(x, y, h) = M^, y) ( ,i\) _ i\ • (3.1) 
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FIG. 1. (i) Cross-section of the vortex lattice on the x-y plane, showing the path along which 
two neighboring lines may cross. The ground state positions of vortices are marked by the small 
circles, (ii) Schematic side-view of the vortex lattice with two lines crossing. The solid lines are in 
the plane of the paper and the dotted lines represent vortices just behind/in front of this plane. 



The subscript, 2, indicates the number of lines allowed to deviate from their positions in 
the ground state triangular lattice. It is now necessary to find an expression for the free 
energy change of the displaced lines compared to the undistorted lattice as a functional of 
a{h), then minimize this functional to find the correct shape of the crossing lines (i.e. the 
configuration of lowest free energy). Looking at ( |2.12D and ( p. 131) one sees that we need to 
expand the differences \ip2\ — iV'ol , and \ip2\ — |^o| , as well as \dip2/dh\ in terms of a. 
The coefficients of these expressions are then integrated over the x-y plane which leads to a 
free energy per unit length along the /i-direction of the form: 
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with the coefficients q equal to the cIq given in Table |, and c^^-* ~ 94.20//^. For example, 
the coefficient c\ was found by calculating: 
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(3.3) 
(3.4) 



Note that all of the integrals that make up these coefficients must be convergent, as the 
denominator of the integrands always grows (for large z) as at least two powers of \z\ higher 
than for the numerator (there is no logarithmically divergent coefficient as the symmetry 
condition ( p.l4| ) is satisfied). The first term in ( p.2|) represents the potential energy cost per 
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FIG. 2. The potential term, J2i H " % fo'^ two lines crossing. 



unit length with the two lines displaced but still straight. The form of this term is shown 
in Fig 1^. The second term is the tilt energy of the two lines. It increases with the distance 
a, which is consistent with an attractive force between the anti-parallel components of the 
vortex segments in the two lines. One could think of the equation ( p.2|) in analogy with a 
Lagrangian density in classical mechanics, of the general form C = T — V for a body moving 
in one dimension, a, and with 'time' equivalent to the /i-direction. The total free energy is 
proportional to the integral over all h of (|3.2|) . Applying the Euler-Lagrange equation for 
stationary values of a functional, and integrating once, one arrives at the equation for the 
form of a{h) that minimizes the free energy change, /2 — a'{df2/da') = const. Substituting 
(p. 21) and applying the boundary conditions gives 
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This can now be integrated up from the point of crossing, a = and (say) h = 0, to give 
the correct form of a{h). As h —^ this will tend to a{h) ~ (4co /c^'^^Y^^y/h (this form at 
small a has been found beforeU). The result is shown in Fig^. The stationary form of a{h) 
can now be substituted into ( p. 2]) and (|2.12| ), and the integral over h performed to find the 
total free energy change for crossing. 

We find / dh f{a{h)} ~ 2 ■ 32, which gives 
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To put this in a simpler form, we use the dimensionless factor ax defined by 



( (3KeiioHkBTyJ2mc\ "" 
"^ = [ ^^ J "^■ 



(3.7) 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



h 

FIG. 3. The form of a{h) that minimizes the free energy subject to the crossing boundary 
condition that a(0) = 0. The expected form as /i ^ is shown. 



Substituting ( p.7| ) into ( |3.6| ) gives 

AFax = ■ 58 kBTlarl^ . 



(3.8) 



All of the following free energies will also be quoted in this form. However, there is at present 
a variety of different units used in the literature. For comparison, we can also write a^ in 
terms of the temperature, the field, and the Ginzburg number Gi. (The Ginzburg number is 
a useful parameter that describes the strength of thermal fluctuations of the superconducting 
order parametertll) . We use the definition of Gi given byllll: 
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The factor (l/Svr)^ enters as we change the units of magnetic energy from c.g.s. to S.I. Using 
standard relations for the coherence length and critical fields we find that: 
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t and h are the dimensionless temperature and field t = T/Tc and h = H/Hc2{0) where 
-f^c2(0) denotes the straight line extrapolation to zero temperature of the Hc2 line when the 
field is applied along the c-axis. 

B. Two Lines Crossing with Relaxation of Nearest Lines 



The result in (|3.8|) is only an upper bound on the actual energy for the two lines to 
cross. This is because we have ignored any relaxation that the remaining lines in the lattice 
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FIG. 4. Cross-sectional view of the lattice when two lines are allowed to cross, and the two 
nearest neighbors move in response. 



may undergo in reaction to the crossing so as to decrease the free energy cost of the defect. 
In order to test how good an estimate this result is, we now investigate the extent that 
allowing the surrounding vortices to move may change the free energy. It seems likely that 
it is the lines closest to the crossing region which will move the most in this defect, and 
therefore affect the energy of the crossing the most. As a second approximation to ( p. 81 ), 
a calculation was performed where the two lines closest to the crossing center, above and 
below, were allowed to move symmetrically in response to the crossing (see Fig ^. The new 
order parameter with this allowed motion of the four lines is: 



ipi{x,y,h) =ipo{x,y) 






(3.11) 



The boundary conditions of this problem are that a(±cxD) = 1/2, b{±oo) = l\/3/2, and 
a(0) = 0. Following the general procedure of the first calculation, the free energy per unit 
length as a function of a{h) and b{h) was found to be: 
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The coefficients q^ , Mj, f j, and Wij are given in Table 0. Again the first term can be viewed 
as a potential term due to the interaction between straight lines. Its form in the a-b plane 
is shown in Fig ^. The other terms describe the tilt energies of the lines, with a fairly 
complex dependence on the positions of the lines. Looking at the form of the potential term 
C4(a,6)=E 
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^"^^ in Fig 1^ we find that the interaction of straight lines in the LLL 
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FIG. 5. The potential term, 6*4(0, 6) = J2i,j ^ij a^^b"^^ , for two lines crossing with the two 
nearest other lines also allowed to move. The dashed line shows the actual path taken in the full 
solution. 



Abrikosov lattice is not at all as we might expect. One would think that parallel vortices 
always have a repulsive force between each other, but instead Fig ^ shows the surprising 
feature that if the two crossing lines are placed at the midpoint, the two nearest neighboring 
lines are actually pulled towards the midpoint! (i.e. the free energy is reduced by moving 
towards the center.) It is impossible to explain this effect by any simple two-body interaction 
between vortex-lines. The minimum value of 6*4(0, b) when a = is 6*4(0, 0.73) ~ 0.54 which 
is considerably less than the value when the two extra lines remain at their original positions, 
6*4(0, v3/2) ~ 0.79. Therefore, allowing nearby lines to move will make a large difference 
to the potential terms for two lines crossing. However, one must also take into account the 
bending terms of all the lines involved in the full solution. Now, to find the configuration 
of a{h) and b{h) that minimizes the free energy cost of two lines crossing, we have to solve 
two coupled Euler-Lagrange equations of the form: 



dxj dh \ dx'i 



0. 



(3.13) 



Substituting ( p. 121) into (|3.13|) leads to two coupled equations that were solved numerically 
with a relaxation method, using a software package. The resulting configuration from the 
solution to (p.l3| ) is only slightly different from the configuration when the nearby lines are 
kept fixed. The two nearest lines move in by ~ 0.03 / at /i = (the actual path is shown in 
Fig|^). The difference in energy is very small, with the same answer to the accuracy quoted 
as (|3.8|) . As we had expected the movement of the nearest lines to those crossing to affect 
the free energy change the most, it seems that the initial approximation of not allowing any 
of the surrounding lines to move is quite good, and we are fairly confident that the result 
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( p.8| ) is close to the exact answer. The reason that the relaxation of surrounding lines has so 
little effect on the crossing energy is that it costs too much free energy to bend these extra 
lines over the short distance along h that crossing takes place, so that the minimum in the 
potential term is never reached. This is a consequence of having a lattice of lines rather 
than pancake vortices. 



C. Search for a Two Line Braid 



The other extension for the two-line problem is not to restrict the lines to moving in one 
direction only but to allow for the possibility that the two lines may 'braid' around each 
other (by braid it is meant that the pair of lines twist around each other by half a rotation 
as we move along the direction of the lines). One of the original motivations to calculating 
the crossing energy was that it was thought to be an energy barrier to two topologically 
distinct configurations. For example, if one considers the two distinct braid configurations 
consisting of either a twist by +7r or by — vr, then to go between the two states one has to 
pass through a configuration with two lines passing through each other. In fact, we show 
that there are no stable two-line braid-defects of the LLL Abrikosov solution (although the 
two-line crossing still remains an important energy barrier between larger entanglements). 

We again allow the two neighboring lines to move symmetrically with respect to each 
other, but throughout the x-y plane. The positions of the two lines are given by {X{h), Y{h)) 
and {—X{h), —Y{h)) respectively (see FigO). The new order parameter is given by 
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FIG. 6. Cross-section of the vortex-lattice with two lines allowed to braid around each other. 

The boundary conditions for a braid are that X(oo) = //2, X(— oo) = — //2, F(±oo) = 0. 
With a similar calculation to the one parameter case, the free energy as a function of X and 
Y can be derived: 
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where Qj are given in Table | (Setting Y{h) = 0, X{h) = a{h) reduces ( p.l5| ) to ( p.2| ) ). The 



J2ij=o cfj'X'^^Y'^^ is the interaction term without bending of the hues, 



first term C2(X,r)=Et=o elf 

shown in Fig 0. We could now go on as before in trying to solve the Euler-Lagrange equations 
to minimize the integral of (|3.15| ) subject to the braid boundary conditions. Inspection of 
Fig 1^ reveals that there is no point in this- there can be no braid solution for this case! This 
is because the form of the potential part of the free energy is always increasing as we increase 
Y away from the midpoint (the midpoint is a saddle-point of C2(X, F)). As the bending 
term will increase with (X^ + Y"^) for any given tilt, it is clear that the free energy cost 
of any two line braid will increase continuously as the path of the braid moves away from 
the midpoint, i.e. there are no stationary configurations for the braid boundary conditions 
other than the direct crossing already discussed. 




, , c]J X'^^Y'^^ , for two lines allowed to braid or cross. A 



FIG. 7. The potential term, C2 = X) 
braid will correspond to going from the minimum at (0.5,0) around to the minimum at (—0.5,0). 
The path with lowest potential will be along the X-axis, which corresponds to a crossing defect. 



One might find it extremely surprising that two straight vortices placed at the same 
point will not repel each other along the Y direction as well as the X direction, even in the 
presence of the surrounding lattice. This seems to be another non-intuitive feature of the 
way vortices interact in the LLL. It is a possibility that the absence of a braid solution is an 
artefact of our approximation of not allowing the surrounding vortices to move. However, 
the result of Section [IIIB| suggests that these surrounding lines should not have a great effect 
on the two-line defect. Also, a similar calculation to that in Section |11IB| was performed 
where the two nearest neighboring lines were allowed to move in response to the motion 
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of the two lines along Y. This still gave the same result that the potential term always 
increases as Y increases. 



D. Three- and Six-Line Defects 

Due to the symmetries of the triangular lattice, there are two other forms of braid defects 
which can be handled in the same way as for the two-line defect. Allowing motion of the 
three lines that form the smallest triangle (see Fig |^), or the six lines that form a hexagon 
surrounding one central line (Fig|^). If we assume that in the lowest energy configurations 
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FIG. 8. Cross-section of the vortex-lattice with three lines allowed to braid around each other. 
The center of the defect is labelled Zf. 
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FIG. 9. Cross-section of the vortex-lattice with six lines allowed to braid around each other. 
The position of the central vortex is z^. 

the lines in these defects move symmetrically about their center, we can write the order 
parameter of the defects as: 
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FIG. 10. Cross-section of the vortex-lattice showing the path taken in a three hne crossing 
defect. 



^3{x,y,h) =iPQ{x,y) 



A{x,y,h) = iJo{x,y) 






{{z-z^) 
The free energy costs per unit length are of the form: 



I') 



(3.16) 
(3.17) 



fn{Rn{h), Onih)} = ^ C^lRn'' COS^^ uO^ sin'^ uOn + C^") R^ 



2n-2 



i,j,k 



dRn 

dh 



+ R, 



2n 



'de^: 

dh 



(3.18) 



with n = 3, 6. The constants cljl and cljl are given in Table 



.(3) 



c^"' ~ 317.9//^ and 

(.(6) r^ 117.7//^^. The conditions for a braid are 6'„(oo) = 0, 6'„(— oo) = 27r/n and -R„(±cxd) 
equal to the ground state value. Again, to minimize the total free energy of given defects, 
two coupled E-L equations must be satisfied. The results of the numerical solutions to the 
E-L equations are: (i) There is no braid solution for the three line defect, other than where 
the three lines meet at the center of the triangle (as in Fig|iy). (ii) The free energy cost of 
the three lines 'crossing' is lower than the cost of two lines crossing, the result being: 



AFsx = ■ 51 keTlarl^ . 
(iii) There is a solution for a six-line braid with an energy cost of: 

AFebr= 2 ■ 19 kBT\aT\t 



(3.19) 



(3.20) 



(iv) There is also a solution for the six-line defect where the lines all meet at a point on 
the center line (see Fig 0), and this crossing defect has a lower free energy cost than the 
six-line braid: 
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FIG. 11. Cross-section of the vortex-lattice showing the path taken in a six hne crossing defect. 



AFex = 1-38 kBT\aT 



(3.21) 



It will seem surprising that we have found the energy for three lines to meet at a point 

to be lower than the energy for two lines. This is due to the potential term when three 

0.54, being lower than the potential term for two lines at their midpoint, 



lines meet, 

.(2) 



.(3) 
-000 



Cqo' — 0.79. Note that both of these values are for when the surrounding lattice remains fixed, 
so it is quite possible that relaxation of the surrounding lattice will reverse this inequality 
in the potentials. As was shown in Section |III B| , allowing the surrounding lines to move 
may decrease the potential at the crossing point significantly. However, we also showed in 



Section PIB| that the full solution including the bending terms in the free energy does not 
allow the surrounding lines to relax much, and so we believe that this result of AF^y. < AF2X 
will hold even when all surrounding lines are taken into account. 



E. Discussion on Two-, Three- and Six-Line Braid Solutions 

To compare the calculations for the above defects, and to understand what is required 
for braid solutions to exist it is useful to cast the free energies (|3.15| , |3.18| ) into the same 
form. We do this by making the transformations, ^ = {X + iY) , ^ = i?3^exp(i36'3), or 
^ = Rq^ exp^iGOe) respectively. This gives a free energy per unit length in terms of the 
complex variable ^: 



uah)} = E 4^[MoniHoy + ^^-^ f i ^^^^\ 



ij=Q 

D„(0 + c?(") 



dh I 



dh J 



d^ 



dh 



(3.22) 
(3.23) 



for n = 2, 3, 6. With this substitution, the two- three- and six-line braids all have the 
same boundary conditions: ^(±cxd) = x„ where x„ is real, and ^ goes around the origin 
of the complex plane once, as h goes from —00 to +00. X2 = {1/2)"^, x^, = (a/3//2)^, 
and Xq = l^. In each case, C, = Xn corresponds to the only minimum of Dn{^) (the ground 
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state configuration). In the analogy with classical dynamics, the solutions to these boundary 
conditions which minimize the integral over h of ( 3.23| ) correspond to the motion of a particle 
of mass m = 2d^'^^ in a two-dimensional potential V{x, y) = —Dn{x + iy). 




FIG. 12. The potential term, DalO = Ei,jd\j'[Re{OY[Im{^)y , where ^ = {X + lYf. A 
two-line braid corresponds to starting and ending at the minimum at (0.25,0) while going around 
the origin once. The path with the lowest potential will correspond to a crossing defect. 



It is clear that a braid solution where the lines in the defect do not meet at their center 
will only exist if there is a 'valley path' in the ^-plane for the function -D„(0 which goes 
around the origin, starting and finishing at ^ = x„. Otherwise any braid configuration will 
be able to continuously lower its total free energy by allowing ^ to move nearer the origin 
as it goes around and the only solution will be a crossing defect (where ^ goes to the origin 
and back along the real axis). By looking at the forms of -D„(0 iii Figs [l^, |13] and |14| one 
can see why there is no braid solution for n = 2, 3, as there are no valley paths. For n = 6 
the presence of a braid solution is explained by the almost circular valley path in D^. 

The presence of a valley path is a minimum requirement for a braid solution to exist. 
It does not guarantee that this braid will have a lower energy cost than the corresponding 
crossing defect, as our results for n = 6 show. Even though the height of the potential term 
Dq{S,) is higher at the origin than at any point in the braid path, the six-line defect can still 
lower its total free energy when ^ goes to the origin and back more than when ^ goes around 
the valley as the bending term in (|3.23|) will be lower when there is less of a change in |^|. 
For the braid boundary conditions, this means that for a defect of given length along h, 
the bending term will be lower when the length of the path in the ^ plane is smaller. Thus 
there is a balance here between the potential term which favours the braid defect around 
the valley, and the bending term which favours the crossing defect. In this case the relative 
terms give the lower total free energy to the crossing defect, although stationary solutions 
exist for both defects with a free energy barrier between the two. 



17 




FIG. 13. The potential term, D^i^) = Ei,j 4f[Re{0Y[Im{0Y , where ^ = i^g^exp (i3^). 
Again, the three-hne braid path with the lowest potential will correspond to a crossing defect. 




^ Re^) 



FIG. 14. The potential term, D^i^) = J^ij 4f [Re{C)Y[Im{^)Y , where ^ = -Re® exp {i6e). Now 
there is a 'valley path' in the potential that corresponds to a braid without the lines meeting at 
the center. 



F. The Twelve-Line Braid Defect 



So far in this paper, the symmetry of the configurations has allowed us to reduce the 
problems to just two coupled equations in two parameters, describing the two-dimensional 
motion of a single 'particle'. In order to consider larger braids than the 2- 3- and 6- line 
defects, one now has to include more degrees of freedom. The simplest problem we can 
think of with more than six-lines involved is the 'twelve line hexagonal defect' (see Fig p!5|) . 
We can still take advantage of the symmetry of this configuration by assuming that it always 
costs least free energy when the six lines on the vertices of the hexagon move symmetrically 
about the hexagon's center, and the six lines on the sides of the hexagon also move sym- 
metrically. This allows us to reduce the number of degrees of freedom in this defect to four 
(corresponding to the two-dimensional motion of two particles in the E-L equations). 



/ ■ 
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FIG. 15. Cross-section of the vortex-lattice with twelve lines allowed to move within a defect. 
The dashed line shows a possible braid path. The dotted line shows the path of a crossing defect. 

In the ground state, the complex positions of the 'edge lines' are given by the roots 



0, and the 'corner lines' by 



0, where zi = 21 and 



of (z — z^) — zi 

Z2 = v^/exp(i7r/6). If the complex positions relative to z^ move from (2:1,^2) to the new 

positions (Ci, C2) then the order parameter of the system can be written: 



i)i2{x, y,h) = ipo{x,y) 



{{z 



Ci')((^ 



C2 



{[Z - ZS - Z,^){{Z - Z.,f ~ Z2^) 



(3.24) 



To find the free energy of twelve-line defects, the same procedure is followed where we 
substitute into ( |2.13| ), then integrate over the x-y plane, to get a free energy per unit length 
as a function of the complex Ci(^) and C2{h)- We omit the details of the calculation, which 
requires a lot of space but results in four coupled non-linear second order equations in the 
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real and imaginary parts of Ci^ and ^2^- This system of four E-L equations was numerically 
solved subject to the braid boundary conditions. The actual path of the braid of lowest free 
energy was close to the hexagonal perimeter, and this solution gave a total free energy cost 
of: 



, 3 



AFubr ^ Q ■ 5 keTlarl^ . (3.25) 

A crossing solution was also looked for, where all twelve lines meet at a point on the 



central line (at z = z^), moving along the dotted line in Fig 15. For the six line case this 



form of defect has a lower energy cost than a braid. For twelve lines crossing the result is 

AFi2x ~ 8 ■ 6 ksTlarl^. (3.26) 



This is higher than the twelve line braid free energy. We expect that all larger braids will 
have lower energies than the corresponding crossing defects. 

IV. CALCULATIONS FOR LARGE DEFECTS 

In this section, the energy cost is calculated for an infinite straight screw dislocation (see 



Fig |T6D, and also for two opposite screw dislocations a large distance apart. The results of 
these are used to give the energy cost for very large braids. The single screw dislocation may 
be a first step towards calculating the energy costs of the screw-edge dislocation loons that 
are thought to be important in describing the dynamics of the Abrikosov crystalt3't3. The 
creation/growth of large braid defects has been proposed as a mechanism for the longitudinal 
resistivity in type-II superconductorglilli3. 

A. Infinite Screw Dislocation 



The order parameter of the infinite screw defect in Fig 16 will depend on only one free 
parameter, s{h): 

, , 1 \ , , xSin7r(z - sVti)) , ^ . 

smvr^; 

where the origin of the x-y plane is now taken to be at the ground state position of one of 
the lines involved in this defect. This order parameter replaces the first order zeros of i\}q at 
z = nl with first order zeros at z = nl + s{h). 

The usual procedure is now followed of integrating the difference of the free energy 
density over the x-y plane. All of the integrands are periodic in x, giving terms which all 
diverge linearly with the system size in the x-direction, L^. Also, the bending term contains 
an integral which tends to a constant for large y, while all other terms are convergent when 
integrated over y. This leads to the result for the free energy per unit length as a function 
of s(/i): 

;,{.„),. f:il.e!".„^.(^)+^c<.. (I): (4.2) 
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FIG. 16. (i) Cross-section of the lattice in the x-y plane showing the path of an infinite straight 
screw defect, (ii) Schematic side-view of the same defect. 



S^) 



with c\"' ~ 0.2883, 4 ^ 2.288 x 10"^ and c^") = vr^d^ol ) ^ 9.8696. To find the general 
form of s{h) that minimizes the integral over h of ( [4. 21 ), for any given system size, we make 
the substitution h ^ h = h/^s with ^g = J Ly/l. This gives a free energy density as a 

function of s{h) with no dependence on Ly and a linear dependence on L^. The resulting 
E-L equation in s{h) is independent of L^ and Ly, so increasing Ly increases the scale of the 
screw defect along the /i-direction, ^s, and the resulting free energy cost (found by integrating 
fs{s{h)} over all h with s{h) the solution to the E-L equation), is proportional to Lx\ L. 



y 



AF, ~ 



0-54- 



L. 



13/2 



y II- 



(4.3) 



That is, the free energy cost increases linearly with the number of lines involved in the 
defect (n oc Lj;) as one would expect, but it also has a divergence as the size of the system 
increases in the direction perpendicular to the defect. This is yet another result that is a 
special feature of the long-range vortex-vortex interactions in the lowest Landau level, and 
one that would not occur with short-range two-body interactions as in the London limit. The 
long range dependence of the bending term in the screw defect has important implications 
to the free energy cost of large braids, which are developed in the next section. 

B. Two Opposite Screws 



A similar calculation which removes this extra divergence of the energy with Ly is of two 
'opposite' screws a large distance Wy = n^/Sl from each other. This defect has the order 
parameter: 



ip2six,y,h) = i;o{x,y) 



sm7r(z - ^^^^i - s{h)) sin 71(2; + ^^i + s{h)) 



iiVs 



sin7r(z 



nl\/3 



i) sin7r(2; + 



nl\/3 



^) 



(4.4) 
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where the origin in {x, y) is now at a vortex half-way between the two opposite screws. For 
large enough n, the potential terms in the free energy cost for a given s{h) become just twice 
the corresponding terms for the single screw. The bending term now contains an integral 
that diverges linearly with n rather than Ly. The actual result is 



f2s{sm = j:2^c'fhin'^ 



1=1 



T 



+ 4nV3^c(^) 




(4.5) 



To solve the resulting E-L equation we make a similar transformation as before: h -^ 

~ 1/2 

h = (2n\/3) • This leads to the same E-L equation as for a single screw, and the free 
energy cost of this 'double screw' is: 



AP 



2s 



1-5- 



Lx, 



p/2 



' y ,2 



(4.6) 



As might be deduced from the result for the single screw, the energy of a large double screw 
is proportional to the square root of the distance between the two screws. It also suggests 
that, at least where the LLL approximation is valid, the energy cost of large braids will not 
simply be proportional to the circumference of the braids (or equivalently the radius of, or 
the number of lines in the braid) as has often been suggested. This is shown in the next 
section. 



C. A Limiting Form For Large Braids 

We now show that as a braid of general shape becomes very large, the free energy cost 
has a simple form depending only on the area enclosed by the braid and the length of the 
perimeter of the defect. We consider a large braid involving n vortices (an example would 
be the braid in Fig|T^, where the lines move from z = Ziio z = Zi^i as h increases between 
±oo, with i = l,n and Zn+i = ^i exp (2i7r). 
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FIG. 17. Example of the shape of a very large braid defect in the x-y plane where the vortices 
move from Zj to Zj+i with increasing h. 

Three assumptions are required to make the necessary approximations: (i) The regions 
over which the braid is straight are much larger than regions where the effect of the 'corners' 
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of the braid on the potential term are important, (ii) All lines in the braid move by the 
same distance, s{h) towards their nearest neighbor, i.e. (iW = Zi + {zi+i — Zi){s{h)/l). (iii) 
The braid is everywhere a large distance from the 'center' of the braid. The center, which is 
where we place the origin, O, is defined by O = (1/n) X]-2j- Therefore this condition means 
that \zi\ ^ /. These assumptions allow us to approximate the potential term in the free 
energy cost using the result for a straight screw (which is proportional to the length of the 
screw) and to calculate the bending term in a simple form as follows. 
The order parameter of this general n-line braid is given by 

I ( u\ I ( J\l=iZ-Zi-{zi+i-Zi){s{h)/l) 
iJnbrix, y, h) = ipo{x, y) — , (4.7) 

llj=l ^ ^i 

To find the bending term we need the partial derivative of the order parameter with respect 
to s: 

d^r.,r , Si=l {-(^.+1 - ^.)(V0 nr=l,^, {Z-Z,- (Z,+1 - Z,){S{K)II))\ 

^^^^° ^:a^^) ^'-'^ 

-^oE^^^^^^^^, (4.9) 

where the second line uses condition (iii). We now consider this sum when the complex 
coordinates z = x ^ iy are far from the perimeter of the braid. That is, when \z — Zj\ 3> 1 
for all j, the sum becomes 



di}n 



br 



ds 




(4.10) 

if z inside braid. {a -i-w 

if z outside braid. 



That is, the partial derivative of the order parameter is periodic within the braid and zero 
outside. Although this result does not hold when z is near the perimeter, when we integrate 
over the x-y plane then for large braids the integral will be dominated by the contribution 
where this does hold. The limiting result for the bending term is: 



(Pr 



P 



dipn 



dh 




^^'(l^ol')!^) , (4.12) 



where A is the area enclosed by the braid. Adding the potential term taken from the infinite 
screw calculation gives a free energy cost for the large braid as a function of s{h) as: 

Ur{sm = j:^<^^^ (^) ^4^''i%]'^ (4-13) 



j=i 



I ' \ I J p \dh 



with P the length of the perimeter of the braid. Making the correct transformation in the 



length scale of the defect, h ^ h = hJPl/AA allows us to use the results from Section [IV A 
to find a total free energy cost for a large braid of: 
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AFnbr ^ 1 • 1^^ kBT\aT\-'. (4.14) 

An important caveat for the large defects in this section is that our results will only be 
good for defects of less than a certain size. As in all of our calculations we have used the 
LLL approximation, we have assumed that we may ignore fluctuations in the microscopic 
flux density. This will lead to good approximations of the energy costs only when the size 
of the defects is less than the distance over which the flux density can change significantly, 
which is the magnetic penetration depth. In the x-y plane, the dimensions of the defect 
must be less than Xat- The extent of the defect in the field direction must be less than Ac. 



V. APPLICATIONS 

A. The nature of the vortex state 

The consequences of finding no stable two-line (or three-line) braids, and of a relatively 
high free energy cost for the six-line braid greatly affects our picture of the vortex state in 
thermal equilibrium at and above the irreversibility line. If stable two-line braids of low free 
energy existed they would be created in large numbers by thermal fluctuations at relatively 
low temperatures, and give rise to a system with vortices twisting around each other in the 
lattice, i.e. an entangled state. However, the absence of such low energy braids means that 
extensive entanglement will not occur over a large region of the H-T phase diagram that 
includes the irreversibility line, as we now show. 

We deduce from our results that the braid defect of lowest free energy is the six-line 
hexagonal braid, with an energy cost given in Eq. ( p.20| ). This is because we have shown 
that smaller braids have no saddle point solutions and the lines can just pull through each 
other while continuously lowering their free energy. We use the result ( p.20 ) to make an 



estimate of how many of these braids will be present in a sample of superconducting material 
of dimensions typical to those available. The stable six-line braid defect extends along the 
direction of the vortices a distance L^br ~ 50 ^c, where ^c is the superconducting coherence 
length along the field direction. If we consider a sample with dimensions along the field axis 
of Lc, then in thermal equilibrium, the average number of braids along any six lines that 
surround a given line will be A^efer = (-^c/-^6fer) exp {—AF^br/kBT). The number of six-line 
braids that any one vortex line is involved in will be 6A^6br- 

If we take an estimate of ax on the melting lineQ as ax ~ —8 we find that the exponential 
gives a factor of exp {—AFghr/kBT) ~ e'^^ ~ 10~^^. Taking an extreme lower limit on the 
coherence length of high-T^ superconductors of ^c ~ 10^ ^°m, we would need a sample size of 
Lc ~ lO^^m for there to be an average of one braid defect on any given line! Alternatively, 
if we consider the number of six-line braids in the whole sample, then at a field of 10 T 
there will be approximately one braid defect in 10^ mm^ at the melting line. (Decreasing 
the external field reduces the number of defects if we remain on the same ax = const line) . 

We now consider at what regimes the braid defects will proliferate. Using the above 
estimates we find that, for a sample with Lc ~ 1mm there will be an average one braid 
defect per line at ax ~ —3. Above this line in the H-T phase diagram, braid defects will be 
plentiful and the vortex state will be strongly entangled. Below this line however, the average 
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line will not be involved in such defects and the vortices can be well identified between the 
top and the bottom of the sample. These arguments make it clear that below the melting line 
in the crystalline phase, and for a large region above in the fiux-liquid phase, the vortices in 
type-II superconductors of sizes currently available will not be in a strongly entangled state, 
when at thermal equilibrium. This is contrary to the picture previously presented inu where 
entanglement has been claimed to play a major role in these systems near the irreversibility 
line. Instead, the lines maintain a correlation between their positions at the top and and 
their positions at the bottom of the sample (for a recent experimental realization, seetJ). 
Note that this does not mean the lines are not moving, either by fiuctuations transverse 
to their length, or their motion as a whole when in the liquid state. It simply means that 
the lines do not twist around each other to form long-lived topologically distinct entangled 
states. 

Of course, our results are only strictly valid in the high field regime where the LLL 
approximation holds. However, if braid defects are of such high energy costs here, it is hard 
to see how the Boltzmann factors will come down as we lower the temperature or field. Also 
our calculations are within the low-temperature crystalline state and it may be questioned 
if we can extend the results to the liquid state. These defects depend only locally on the 
surrounding lattice, and we expect the liquid state to retain crystalline order over such 
length scales for a large region above melting, in which case our results may still apply. It is 
also important to note that the absence of entanglement is a size dependent phenomenon. 
In a big enough sample, entanglements will arise in principle, but, at least near the melting 
line, they will be completely absent in practice. The above arguments are only for a system 
in thermal equilibrium, and if the system has been rapidly cooled from a high temperature 
to near the melting line it is still possible for entanglements to be present for some time. 

The absence of any long lived braid defects within the system of vortices can be tested 
experimentally. Neutron scattering, or possibly muon spin resonance, could be used to inves- 
tigate the magnetic fields associated with the braids which are transverse to the externally 
applied field. There should be a locally qualitative difference between the resulting spectra 
due to the fluctuations in the transverse displacement of a vortex along its length, and from 
the local transverse fleld of a stable braid defect. 



B. Longitudinal Resistivity due to Braid Defects 

Experiments have been performed where a current is applied between terminals on the 
top and bottom faces of a superconductor, with the external magnetic fleld applied parallel to 
this currenttlnl3, allowing measurement of the longitudinal resistivity. It has been known for 
some time that the spontaneous production of flnite vortex loops at non-zero temperatures 
can be a source of non-linear resistivity in the Meissner-phase of any superconductorEO. This 
idea has been extended to the problem of the resistivity along the fleld direction in the mixed 
state of a type-II superconductorEJ. If one projects the vortex lines onto the x-y plane, then 
entanglements of the vortices may be identifled as 'planar loops' perpendicular to the fleld 
direction. For any current, Jz, flowing in this direction, these planar loops will have an 
interaction with the current analogous to the vortex loops in the Meissner-phase. The braid 
defects that we have calculated are precisely the sort of defect that project to a planar loop. 

It is important to distinguish between two different limiting cases. The case we consider 
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is where the vortex-lattice contains 'weak entanglement' only, i.e. entanglements produced 
by the presence of the current itself. This should be distinguished from 'strong entanglement' 
where planar loop type entanglements will exist at equilibrium, on all length scales, which 
gives rise to a linear longitudinal resistivityll3. However, we have shown that in a large region 
of the phase diagram, there are very few entanglements present in the size of samples used 
in experiments, so that the systems are always in the weak entanglement limit. The fact 
that the resistivity has been observed to become non-linear below a line much higher than 
the irreversibility line in YBCC)ll3 supports our claims in Section [VA| concerning the absence 
of entanglements. 

Following the arguments of Feigelman et. al.Ej, but using our results for the free energy 
costs of the relevant defects, we find that in the presence of a current density J^, a large 
braid has a total free energy cost of: 

AF(A, J,) = eA^t - 0o^.A (5.1) 

with eA ^ 2.ll-'^/^kBT\aTf'^ = 2.3 {B/cpof^'^kBTlaTf'^. Therefore there will be a crit- 
ical area, Ac{Jz), above which the braid will be driven to larger growth by the magnetic 
interaction with the current: 

This will cause dissipation, which will be proportional to the probability of nucleating large 
enough braids by thermal fluctuations. The nucleation of large enough braids requires 
jumping a free energy barrier, which will lead to a longitudinal resistivity of the form: 

Pzz = -^^e y'^J , (5.3) 

with Jt — 1.43 (-B/0O )kBTaT^. This form is quite different to that given previouslyt^, 
where it was assumed that the free energy cost of planar loops had a linear dependence on 
their radius. 

Unfortunately, the applicability of the above results to any experiments that measure 
this longitudinal resistivity is rather doubtful. This is because at regions of the H-T phase 
diagram of interest, and for all reasonable current densities, the critical radius of a braid that 
needs to be created is far larger than the magnetic penetration depth Aab, (i.e. the length 
scale in the ah plane over which the magnetic field may vary significantly). Therefore, the 
important braid defects will be outside the scale where the LLL approximation may hold. 



C. Crossing Defects as Energy Barriers to Disentanglement 

The energy barrier t/x to the crossing of two vortices, or to a process of cutting and 
reconnection, is the energy cost of the configuration with the highest free energy that must 
exist during these processes. This configuration is assumed to be where the two lines meet 
each other at a point, as in the crossing defects that we have calculated. Therefore, our result 
for AF2X for the free energy cost of a two-line crossing defect maybe used as an estimate of 
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Ux ■ The value of Ux has been shown to be an important parameter of the entangled vortex 
stateffl. The relaxation time of the entangled vortices grows exponentially with Ux/ksT. 
A different form for the longitudinal resistivity to ( |5.3| ) that applies when there is strong 
entanglementtJ (i.e. braids exist on all length scales) is proportional to exp(— [/x/ZcbT) (as 
crossing defects will be the basic energy barrier to the growth of braids). It is doubtful 
where this applies though, as we have shown above that the vortex state is not entangled 
over a large region of the H-T phase diagram. 

Estimates of Ux have also been usedH to explain the transition from local to non-local 
resistivity seen in the DC flux transformer experiments on YBCOta'EI. In this interpretation, 
the average length over which crossing defects take place is associated with a correlation 
length of vortex 'identity H When this length scale becomes smaller than the system size, 
the vortices will be uncorrelated between the top and bottom of the sample, and non-local 
effects will be of less importance in the resistivity. This argument does not depend on the 
entanglement of the vortices. 



VI. SUMMARY 

To summarize, we have calculated the free energy costs of an assortment of topological de- 
fects of the Abrikosov vortex-lattice that is the mean-field solution for the Ginzburg-Landau 
free energy functional of a perfect bulk Type-II superconductor. All of our calculations are 
made within the lowest Landau level, which is a good approximation at high fields just below 
the Hc2 line. Although we have made a large approximation in not allowing the surrounding 
lattice to move in response to these defects, we believe that our answers are still very close 
to the full solutions with relaxation included. This is a consequence of the elastic line nature 
of the lattice, whereby the vortices' resistance to tilting holds them near to their ground 
state positions. 

The free energy costs calculated are for crossing defects of 2, 3, 6, and 12 lines; braid 
defects of 6 and 12 lines (see Table [TVD ; infinite screw defects and a form for very large 
braids. Within the LLL, one cannot truly view the vortex system in terms of discrete lines 
with simple two-body interactions between them (as is the case in the opposite extremes 
of the London limit). This non-trivial nature of vortex- vortex interactions has led to a few 
unexpected results of our calculations. Firstly there were no braid solutions for two- or three- 
line defects - if any braid configuration of two/three lines is imposed it may continuously 
lower its energy by the lines meeting each other at a point (and thus forming a crossing 
defect)! Secondly, we have found a lower free energy cost for three lines crossing than for a 
two-line crossing. One assumes that this is a consequence of the lattice symmetry, together 
with the effect stated above that the surrounding lines in the lattice do not relax much in 
response to a localized defect. Thirdly, it costs a lower free energy for a six-line crossing 
defect than a six-line braid, even though the potential terms are higher for the crossing 
defect (the tilt terms decrease when the six lines move nearer each other). This odd result 
was not carried through to the twelve-line defects, and we believe that all larger braids cost 
less energy than the associated crossings. Finally, we found that for large screw defects, the 
free energy cost diverges with the system size perpendicular to the defect. This effect (a 
long range dependence in the tilt energy) leads to a result for the free energy cost of large 
braids which depends on the area enclosed by the braid in a new way. 
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We think that our results will be useful in interpreting experiments on Type-II supercon- 
ductors where the dynamics of the vortices play an important role, at least in the regimes 
of high field, and where any layered structure of the superconductors may be ignored (i.e. 
as long as the coherence length ^^ is greater than the inter-layer thickness). Some particular 
areas of application were outlined in Section |V], but there are still some remaining questions 
yet to be cleared up. For instance, what if any are the consequences of a lower energy 
cost of crossing for three-lines than two-lines? We have not yet been able to use the en- 
ergy costs calculated to make quantitative predictions on transport properties such as the 
resistivity. However, we have shown that the equilibrium vortex state will not be entangled 
over a significant region above and below the irreversibility line of a high-T^ superconductor. 
This prediction is open to experimental test, possibly with a magnetic probe like neutron 
scattering. 
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APPENDIX A: PRODUCT FORM OF THE JACOBI-THETA FUNCTION 

In this appendix, it is shown explicitly that f{z) = ^^{ttz/I^ttt/I) has the form of 
This then justifies our procedure in forming a new function with displaced zeros, but 
which is still in the LLL subspace, as is done in ( f^.ll| ). Starting from the standard product 
representation of a Jacob i-theta functionO, and for clarity letting / = 1, we have 

oo 

^3(vr2, Trr) = e n (l - e^'^^"""")) (l - e'^^'^'^'-^) , (Al) 

ra=0 

where 9 = H^o (-'- ~ e^'^*"'^) = const, and Zn = —{n + |)r — i. (For the triangular lattice, 
r = 1/2 + i\pij2 and O ^ 1.0043.) This is the form that was used to find i\}q when 
numerically calculating the coefficients of Tables |-|I|. Now, the two products in (|A1|) can 
be rewritten as (1 — exp [±27rz(2; =f 2;„)]) = =i=2zexp [±7ri(z =F z„)] sin7r(z =]= 2;„). Substituting 
this expression into (^l|) gives 



^3(7r2,7rr) = J| 4e"^*''^" sin7r(2 - z„) sin7r(2 + 2„). (A2) 

n=0 

The prefactor, H^o^^^P (~2mz„), can be seen to have the correct form to ensure conver- 
gence of the product, because as n — > cxo, for finite z, we get: sin7r(2 — z^ sin 71(2; + z„) ~ 
(1/4) exp (2i7rz„). Therefore for large n, the terms inside the product tend to unity (as is 
necessary for an infinite product). Now we can use the well known product representation of 
sinx = xO^i (1 ~ (x^/n^TT^)). Putting the sine terms in ( [A^ ) in to this form finally reveals 
the product form we are looking for: 

^3(vrz,7rr) = 7r^0n4e-^^-'(z + z„)(^-z„) fl f 1 " ^^^^) f 1 - ^^^^V (A3) 

n=0 ij=l \ ^ / \ J / 
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Each of the terms in the second product may be expanded as 1 — {z ± Zn) j 3^ = 
(1/i^) (i + (yZ ± Zn)) (j — {z ± Zn)). This is proportional to the product of two terms of 
the form {z — Zj), where the roots, Zj, are given by 

Zj = TZn ± j 

= ± (^^ + nr] ± j. (A4) 



Therefore, we have shown that f{z) = i)^ {'n'z/l, tit/I) has the general product form of 
and the zeros are positioned periodically in the complex plane, with periods r and 1. 
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TABLE I. Coefficients for two lines allowed to move, c, 



(2) 



1 





0.7905 


-5.201 


10.0395 


-32.7365 


100.8234 




5.201 


5.1655 


-32.7365 


403.2937 




J2) 
Cj2 


10.0395 


32.7365 


604.9406 






Cj3 


32.7365 


403.2937 










100.8234 











TABLE IL Coefficients for two lines crossing with two lines relaxing, q- , Ui, Vi, Wij 
12 3 





2.0304 


6.2416 


8.6774 


5.9207 


3.3321 




-6.2416 


-19.9304 


-27.4202 


-33.9090 


-7.5801 


J4) 
Cj2 


8.6774 


27.4202 


40.8076 


39.0456 


48.8558 


Ci3 


-5.9207 


-33.9090 


-39.0456 


-103.4540 


-26.7255 


J4) 
Cj4 


3.3321 


7.5801 


48.8558 


26.7255 


274.8685 


Ui 


28.7674 


-13.9138 


28.7674 






t'j 


28.7674 


13.9138 


134.871 






WiO 


57.5348 


-13.9138 








Wil 


13.9138 


269.742 









TABLE III. Coefficients for three and six lines allowed to move, c-A, clA 


ijk 


000 200 400 110 220 202 


310 




0.5443 -23.7064 293.7881 -2.0347 11.6525 11.5790 
1.40174 -2.8908 1.32786 -0.18099 0.37540 0.38420 


-15.3239 
-0.03321 


TABLE IV. Summary of the Free Energy costs of Crossing and Braid Defects 




i 2 3 6 


12 


AFix/(A;BT|aT|2) 0.58 0.51 1.38 
AF,br/{kBT\aT\^) 2.19 


8.6 
6.5 
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